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Recently it was realized ^ that linear optics and photo-detectors with feedback can be used for 
theoretically efficient quantum information processing. The first of three steps toward efficient linear 
optics quantum computation (eLOQC) was to design a simple non-deterministic gate, which upon 
post-selection based on a measurement result implements a non-linear phase shift on one mode. Here 
a computational strategy is given for finding non-deterministic gates for bosonic qubits with helper 
photons. A more efficient conditional sign flip gate is obtained. 



I. INTRODUCTION 

Now that we know that linear optics and photo-detectors are sufficient for quantum information pro- 
cessing [|l|], it is necessary to investigate how the necessary schemes can be realized more efficiently. One 
promising direction is to use superpositions of squeezed or coherent states for encoding qubits |^ ^]. In this 
note, it is shown how the non-deterministic gates at the foundation of the constructions in [|I|] can be found 
and improved. Other relevant work in this direction includes [Q, ||, |6|, ^], where networks suitable for ex- 
perimental realization are given. The focus of this paper is on what can be done in principle without giving 
experimentally accessible layouts. To that end, a systematic method is given for finding non-deterministic 
gates based on a combination of algebraic solution finding, exploitation of known symmetries, and nu- 
merical optimization. By using the method, a conditional sign flip for bosonic qubits that succeeds with 
probability 1/13.5 using two helper photons is found. This improves the one in which succeeds with 
probability 1/16. What is the optimum probability of success for any number of helper photons? A charac- 
terization of the achievable states without post-selection implies that the probability cannot be one, a result 
related to known bounds on Bell-measurements ^ 

II. PRELIMINARIES 

The physical system of interest consists of bosonic modes, each of whose state space is spanned by 
the number states |0), |2), If more than one mode is used, they are distinguished by labels. For 
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example, \k)^ is the state with k photons in the mode labeled r. The hermitian transpose of this state is 
denoted by '^{k\. The vacuum state for a set of modes has each mode in the state |0) and is denoted by |0). 
The anniliation operator for mode r is written as a^''^ and the creation operator as c^''^ = (a'^'^^)^ Recall 
that c('')|m) = + l\m + 1). Labels are omitted when no ambiguity results. Hamiltonians that are at 
most quadratic in creation and annihilation operators generate the group of linear optics transformations. 
Among these, the ones that preserve the particle number are called passive linear. Every passive linear 
optics transformation can be achieved by a combination of beam splitters and phase shifters. If U is passive 
linear, then Uc^^^\0) = UsrC^^^IO), where Ugr defines a unitary matrix u. Conversely, for every unitary 



matrix u, there is a con^esponding passive linear- optics transformation |1C]. For the remainder of this note, 
all linear optics transformations are assumed to be passive. 



m. CONDITIONAL PHASE SHIFTS 

A conditional phase shift by 9 on two modes is the map CSg : \ab) — > e*("^)^|a6) for < a, 6 < 1. 
These phase shifts can be used to implement conditional phases on two bosonic qubits. A bosonic qubit 
Qr,s is defined by identifying logical |o) with |01)^^ and logical |i) with |10)^^. The modes r and s can be 
two distinct spatial modes or the two polarizations of one spatial mode. To realize the conditional sign 
flip between Qi 2 and Q3,4, apply CSi8o° to modes 1 and 3. The bosonic qubit controUed-not can then be 
implemented using conditional sign flips and single qubit rotations, which are realizable with beamsplitters. 

In conditional sign flips were implemented indirectly using a non-deterministic realization of the 
map 

NS : a|0) + + j\2) ^ a|0) + - 7|2) (1) 

that succeeds with probability 1/4. This realization requires one helper photon and two ancilla modes. 
The goal is to implement CS51 more efficiently directly using two helper photons. One helper photon can 
be shown to be insufficient by means of the same algebraic method about to be used. Let modes 1 and 2 
contain the state to which CSq is to be applied. The basic scheme is to start with two ancilla modes 3 and 
4 initialized with one photon each, apply a linear optics transformation to modes 1,2, ... ,k with k > A, 
measure all but the first two of these modes and accept only a predetermined outcome, say where one photon 
is detected in each of modes 3 and 4 and none in the added modes. Let u be the unitary matrix associated 
with the lineal- optics transformation, with Usr the entries of u. The post-selected final state is determined 
completely by the 4x4 upper left submatrix F of u with entries Vrs = Ugr for s, r < 4. 

It is necessary to consider the effects of the scheme on the initial states 1 00) , |01) ,|10) ,|11) .Since 
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photon number is conserved, we have, without renormalization: 



1 00) ^2 ^ aoooo|00)^2 (2) 

101) ^2 ^ "0101101)^2 +"oiio|10\2 
110)^2 ^ aioiol 10)^2 +aiooi 101)^2 

\l\ ^ aiiii|ll\2 + aii2o|20)^2 + «ii02|02\2 

To be successful, the amplitudes have to satisfy 

"0110 = aiooi = "1120 = "1102 = (6) 
"1010 = "0101 = "0000 (7) 
"nil = e*^"oooo (8) 

The amplitudes are polynomials of the coefficients of V. For example aoooo = ^33^^44 + ^34^^43- More 
generally, define ps = VrsC^^^ ■ If the initial state in mode s has dg photons, then the output state is given 
by P = Pg". Thus, P is a polynomial of the c^''^. If (3 is the coefficient of the monomial (c^*))"^' in 



P, then the output amplitude for having nit photons in mode t is given by ^711* mtW- This shows that the 
amplitudes Oabcd are polynomials of the Vrs- 

The first step for constructing CSg is to solve Eqs. |6|-|8|, which are polynomial identites in the Vrs- Before 
showing how to reduce the difficulty of doing that, let us see how to proceed from there. Since there are 
16 free complex variables, the solution will have a number of remaining free variables that must be chosen 
to optimize the probability of success ( | aoooo P) and to satisfy one more constraint: The solution is an 
(explicit) matrix V that needs to be extended to a unitary matrix u. This is possible if and only if the 
maximum singular value (that is the square root of the maximum eigenvalue of V'^V) is at most one. The 
extension is not unique. One can set the first four columns of u to the matrix with orthonormal columns 

^=1^1, (9) 
\(I-yty)i/2 J 

and then complete the last four columns with any orthonormal basis of the orthogonal complement of the 
space spanned by the columns of X. The maximum singular value constraint is needed to be able to compute 
the square root in the expression for X. If some of the singular values of V are equal to one, then fewer 
than four additional columns and rows can be used 

The singular- value constraint cannot be easily achieved using algebraic methods. In principle, one can 
reparametrize the matrix V to guarantee the constraint, for example by using the polar decomposition and 
an Euler angle representation of unitary matrices. In the case where CSg is to be applied to the "left" 
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modes of a pair of bosonic qubits, the singular value constraint can be removed by exploiting a rescaling 
symmetry. Now there are two additional modes to complete the bosonic qubits. The total number of photons 
is always four. Let F be a matrix whose coefficients satisfy the identities for the Oabcd- Let A = X{V) be 
the maximum singular value of V and consider the matrix 



where the first two indeces are associated with the two other ("right") qubit modes. Ve has maximum 
singular value 1 and can be extended to a unitary Ug as before. The claim is that if the resulting optics 
operation is applied to the pair of bosonic qubits with the same post-selection procedure, it has the intended 
effect with probability To see that this is true, first observe that, V = XVe satisfies the polynomial 

equations obtained by requiring that the operation works correctly for the pair of bosonic qubits. The 
amplitudes (as in Eqs. ^|-^ that occur in these equations are polynomials which ai^e either homogenous 
linear in the coefficients of a given column of V' or independent of them. This is because the input states 
have at most one photon in each mode. Because each input state under consideration has exactly four 
photons, the amplitudes are all homogenous of degree four in the coefficients of V'. This implies that 
multiplying V' by 5 scales the amplitudes by 5^. Since the equations to be satisfied are homogenous linear 
in the amplitudes, every scalar multiple of the matrix also satisfies the equations. 

With the observation of the previous paragraph, instead of trying to satisfy the singular value constraint, 
one can recalculate the probability of success by dividing V's probability of success by and optimize it. 
Note that this works even if A < 1. In computer experiments using naive optimization methods (see below), 
this usually led to solutions V with A = 1 for 6* = 180° and 9 = 90°. 

To simplify solving the equations for the Uabcd^ one can use scaling symmetries to standardize V. Since 
each of the Uatcd is homogenous in the variables associated with any one row or column of V, the equations 
of the form Oabcd = are satisfied for any rescaling of a row or column. The non-zero aabcd ai'e homogenous 
of degree one in the third and fourth column (because of the presence of the helper photons at the input in 
modes 3 and 4) and in the third and fourth row (because of the post-conditioning on detecting exactly one 
photon in each of modes 3 and 4). Because Eqs. ^|-|8| are homogenous linear- in the amplitudes, rescaling 
these rows or columns preserves the identities. The non-zero aabcd satisfy that they are of equal degree and 
homogenous in the first column and (separately) in the first row. This is due to the fact that when a photon 
is present in mode 1 at the input, this is designed to be the case at the output too. Thus multiplying the first 
column by 6 and the first row by 1/5 does not change the values. Similarly, this rescaling can be used on 
the second column and row. 
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The scaling rules of the previous paragraph can be used to introduce unconstrained scaling variables and 
standardize the entries of V. For example, one can take 1^13 = ^24 = U33 = f44 = ^43 = 1. Note that this 
choice implies that solutions where any one of these variables is are not easily found. It may therefore be 
necessary to try solving with some of the variables set to 0. For example, the CS gate of [|l|] after translation 
into the form used here satisfies ^43 = 0. I did not find any solutions satisfying this constraint with better 
probability of success. 

Mathematica was used to solve the equations (any other computer algebra system would do equally 
well). The strategy was to solve linear equations first and then to simplify expressions. Some Mathematica 
notes are included verbatim in App. ^ and include formulas for the solution found. The solution could be 
expressed in terms of the remaining variables of the last two columns of V and one additional variable. 
After some experimentation, it seemed that in all the best solutions, vu = V22- This was exploited in the 
final version of the optimization procedure, implemented in Matlab (the programs are available by request). 
Briefly, the function to be optimized takes as input the remaining free complex variables (^14, ^23, V34, ^i), 
and a non-redundant subset of the scaling variables. To avoid infinities, one can provide the logarithms 
of the scaling variables as inputs. The scales can be taken to be real since phases have no effect on the 
probability of success. The function can then be optimized using random starting points. With the opti- 
mization procedures provided by Matlab, it was found useful to randomly perturb the point returned and 
repeat until the solution no longer changes significantly. This procedure routinely finds the same optimum. 
For 9 = 180°, it was possible to guess the algebraic numbers that it converged to. Here is a version of the 
matrix found, which turns out to be unitary: 

-1/3 -\/2/3 V2 /3 2/3 ^ 

V2/3 -1/3 -2/3 V2/3 

^ ^ , / (11) 

V3+V6/3 V3 - \/6 /3 -•^(3 + ^/6)/2/3 ^1/6-1/(3^/6) 

VS-VG/S -\/3 + V6/3 -^1/6-1/(3^6) --^(3^7^/3 y 

The probability of success is 2/27. The matrix can be systematically decomposed into elementary beam 
splitter and phase shift operators [[lC[]. An optical network realizing it is shown in Fig. [l[ The implementa- 
tion uses fewer elements (4 beam splitters, 4 modes, 2 photon counters, probability of success 2/27) than 
the solution in [|l|] (6 beam splitters, 6 modes, 2 photon counters, 2 photo-detectors, probability of success 
1/16). As before, the counters must be able to distinguish between zero, one, or more than one photons. 



180° 




FIG. 1 ; Optical network realizing CSi8o° ■ The notation is as explained in [|l]] . 



90° 
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A matrix that can be extended to obtain CSgo" by post-selection was also obtained: 

/-0.3202 + 0.0418i -0.2520 - 0.3226i 0.2883 -0.1292 - 0.722lA 

-0.2520 - 0.3226i -0.3202 + 0.0418i -0.1292 - 0.7221z 0.2883 

-0.3216 + 0.7210i -0.1711 - 0.1725i 0.2469 0.3322 -h 0.3285i 

\^-0. 1711 - 0.1725i -0.3216 -h0.7210i 0.3322 -F0.3285i 0.2469 / 

The probability of success for this solution is 1/19.37. A "nice" beam splitter decomposition of this matrix 
was not found. Partly this is due to the fact that because only two of the singular values are (close to) one, 
at least two extra modes must be added for the unitary completion. The simplest method of decomposing a 
6x6 unitary matrix normally requires 15 beam splitters. 

It is an open problem to determine whether the above solutions are indeed optimal as is suggested by the 
results of the numerical experiments. 



IV. BOUNDS ON CONDITIONAL PHASE SHIFTS? 



To obtain bounds on the probability of success of a phase-shift gate implemented with helper photons, 
one can attempt to characterize the states obtained in the output modes after tracing out the helper modes. 
There is some choice of the initial state of the modes that the gate is applied to. Assume that this is a state 
obtained by applying linear optics to prepared single photons. In this case, the final state after a linear optics 
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transformation is given by 

IV'/) = U^W^''^ + ... + «fc,nc('^))|0). (13) 

fc=l 

The goal is to show that after tracing out modes m' + 1, . . . , m, the state in the remaining modes is a mixture 
of states of the form 

n 

HiPkO + PklC^'^ +■■■+ Pkm'C^'^'^m. (14) 
k=l 

In fact, this is the case if the final state before tracing out is also of this form. To be explicit, add to the 
factors in the expression for \ipf) any constant terms ako so that 

n 

l^/> = Uio^ko + afcic(') + . . . + afc^c("^))|0). (15) 

k=l 

First trace out mode m. Given an overcomplete set of states |7)^, the state of modes 1 ... m is a mixture of 
the (unnormalized) states '^(7||V'/)- Choosing as the set of states the coherent states and using the fact that 
for these states '^(7|c('^) = "^(717, the mixture consists of states of the form 

71 

JJ(afco + 7afem + afcicW + --- + afc(m-i)c("'"^^)|0)^ . (16) 

k=i 

Iterating this procedure proves the desired result. 

Consider the conditional sign-flip gate. With this gate and using a few beam splitters, one can map the 
state |1100) to the state -^{\nOO) + |0011)) = -^^(c^^^c^^) + c^^^c^'*)), a well-known entangled photon 
state. By the above, before post-selection on a measurement of the other modes and with n helper photons, 
the state can be written as a mixture of products of linear expressions in the creation operators. To obtain 
a bound on the probability of success, it suffices to obtain a bound for the overlap of (normalized) such 
states with the Bell state. Because the normalized overlap of (c^-^^ + c('^))(c(^) + c^^^) with the Bell state is 
1 / \/2, the bound on the probability of success thus obtained can be no smaller than 1/2. It is clear that the 
probability of success cannot be made equal to one: The polynomial xy + uv associated with the creation 
operators in the Bell state cannot be factored. 

A problem suggested by the above is: 
Problem. What is the maximum probability of success for implementing CSg using linear optics with at 
most k independently prepared helper photons and post-selection from photon counters without feedback? 

It was shown that for 9 = 180°, a probability of success of one is not possible, but for k > 2, 1/13.5 
can be realized. A variant of the problem asks the same question for conditional sign shifts of two bosonic 
qubits (a four mode operation). Other directions for investigation are to determine what improvements are 
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possible if active linear optics operations can be used, or if initial states such as prepared entangled photon 
pairs ||^ or photon number states like |2) are available. 

Acknowledgements. This work was partially supported by the NSA and the DOE (contract W-7405-ENG- 
36). 
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APPENDIX A: MATHEMATICA NOTES FOR SOLVING EQ. @-| 

The following gives a sequence of steps using Mathematica for solving Eq. |^-|8[ 

{* The columns of V as polynomials of the creation operators 
* cl,c2,c3,c4 with the scaling simplifications. *) 
pi = {Vll*cl+v21*c2+v31*c3+v41*c4) ; 
p2 = (Vl2*cl+v22*c2+v32*c3+v42*c4) ; 
p3 = {cl+v23*c2 + c3 + c4) ; 
p4 = {vl4*cl + c2+v34*c3 + c4) ; 



{* The amplitudes that occur in the equations to be solved are: *) 



aOOOO := 


Coef f icient [p3*p4, c3*c4] 






alOlO := 


Coef f icient [pl*p3*p4, c3* 


c4 


* c 1 ] ; 


alOOl := 


Coef f icient [pl*p3*p4, c3* 


c4 


*c2] ; 


aOllO := 


Coef f icient [p2*p3*p4, c3* 


c4 


*cl] ; 


aOlOl := 


Coef f icient [p2*p3*p4, c3* 


c4 


*c2] ; 


aim : = 


Coefficient [pl*p2*p3*p4. 


cl 


*c2*c3*c4] ; 


all20 := 


Coefficient [pl*p2*p3*p4. 


cl 


"2*c3*c4] ; 


all02 := 


Coefficient [pl*p2*p3*p4. 


c2 


"2*c3*c4] ; 


{* First 


solve al001==0. 







* alOOl is linear in the coefficients of pi. *) 
plx34 = Coef f icient [alOOl, v21]; 

plx24 = Coef f icient [alOOl, v31]; 
plx23 = Coef f icient [alOOl, v41]; 

{* Parametrize the linear solutions by hand, introducing 112,113,114 
rll = { v21->112*plx24 + 113*plx23, 

v31->112* {-plx34) + 114*plx23, 

v41->113* (-plx34) + 114* (-plx24) } ; 
pi = pi/ . rll ; 

{* Similarly for a0110==0. *) 
plx34 = Coef f icient [aOllO, vl2]; 
plx24 = Coef f icient [aOllO, v32]; 
plx23 = Coef f icient [aOllO, v42]; 
rl2 = { vl2->121*plx24 + 123*plx23, 

v32->121* {-plx34) + 124*plx23, 

v42->123* (-plx34) + 124* (-plx24 ) } ; 
p2 = p2/ . rl2 ; 

{* Next solve al010==a0101==aOOOO : 

* Again this leads to linear equations in vll and v22 respectively. 



xl = Coef f icient [alOlO, vll] ; 
aplOlO = Simplify [al010-xl*vll] ; 
rl3 = {vll-> (aOOOO-aplOlO) /xl } ; 
pi = pi/ . rl3; 

x2 = Coef f icient [aOlOl, v22] ; 
apOlOl = Simplify [a0101-x2*v22] ; 
rl4 = {v22-> (aOOOO-apOlOl) /x2 } ; 
p2 = p2/ . rl4 ; 

(* Now solve for all20==0 and all02==0. 

* This leads to quadratic equations in 112 and 121. 

* Explicitly: *) 

1111 = Coef ficient [FullSimplify [all20] , 112*121] ; 
tlll20 = all20-llll*112*121; 
1101 = Coef ficient [tlll20, 121] ; 
1110 = Coef ficient [tlll20, 112] ; 

1100 = FullSimplify [tlll20-1101*121-1110*112] ; 

Imll = Coef ficient [FullSimplify [all02] , 112*121] ; 
tlll02 = all02-lmll*112*121; 
ImOl = Coef ficient [tlll02, 121] ; 
ImlO = Coef ficient [tlll02, 112] ; 

ImOO = FullSimplify [tlll02-lm01*121-lml0*112] ; 

xysols = Solve [x*y*llll + x*1110 + y*1101 + 1100 == 
x*y*lmll + x*lmlO + y*lm01 + ImOO == 0, {x,y}]; 
xyl = FullSimplify [xysols [ [1] ]] ; 

pi = (pi/ . {112->x, 121->y} ) / .xyl; 
p2 = (p2/ . {112->x, 121->y} ) / .xyl; 
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{* 

* One can now simplify the expressions by removing 

* redundant variables introduced earlier, pi's coefficients 

* are a function of Coefficient [pi , c4 ] and similarly for 

* p2. *) 

Isimrule = { 113-> (11-114* {l+v23) )/ {l+v34 ) , 

123-> (12-12 4* (1+V14) ) / (l+v34) } ; 
pi = pi/ . Isimrule; 
p2 = p2/ . Isimrule; 

(* Check the identities by evaluating: 

* Answers included after the expression: 
FullSimplify [aOOOO] //InputForm 

1 + v34 

FullSimplify [aOlOl] //InputForm 
1 + v34 

FullSimplify [alOlO] //InputForm 
1 + v34 

FullSimplify [aOllO ] //InputForm 


FullSimplify [alOOl] //InputForm 


FullSimplify [all02] //InputForm 


FullSimplify [all20] //InputForm 


FullSimplify [allll ] //InputForm 

(-1 - 4*ll*12*v23"2*vl4~2* (-1 + v34)*v34 + 

v34*(-l - 4*11*12*(-1 + v34) + v34 + v34"2 + 2*v23*(l + v34)"2) + 
vl4*(-2*(l + v34)"2 + v23*(-l + v34)*(l + v34*(2 + 8*11*12 + v34))))/ 
((1 + v34)*(-l - (2 + v23)*vl4 + v34 + v23* (2 + vl4)*v34)) 
*) 



(* Remaining identity: allll = ph*aOOOO, where ph is the desired phase. 
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* Note that allll is now a function of 11*12, so solve for that. *) 
llrls = Solve [FullSimplify [allll == ph*aOOOO], {11}]; 

llrl = FullSimplify [llrls [ [1] ]] ; 
lltl2 = 11*12/ .llrl; 

{* 

* Checking shows that pl-cl and p2-c2 are multiples 

* of 11 and 12 respectively. Exploit that to express 

* the coefficients of pi and p2 : *) 

dplcl = FullSimplify [ (Coefficient [pi, cl] -1) /II] ; 
dplc2 = FullSimplify [Coefficient [pi, c2] /II] ; 
dplc3 = FullSimplify [Coefficient [pi, c3] /II] ; 
dplc4 = FullSimplify [Coefficient [pi, c4] /II] ; 

dp2cl = FullSimplify [Coefficient [p2, cl] /12] ; 
dp2c2 = FullSimplify [ (Coefficient [p2, c2] -1) /12] ; 
dp2c3 = FullSimplify [Coefficient [p2, c3] /12] ; 
dp2c4 = FullSimplify [Coefficient [p2, c4] /12] ; 

(* The transpose of V is now given by: 
*) 

vmat = { {dplcl, dplc2, dplc3, dplc4 } *ll+{ 1, 0, 0, } , 

{dp2cl, dp2c2, dp2c3, dp2c4 }*12+{0, 1,0,0}, 

{l,v23, 1, 1}, 

{vl4, 1, v34, 1} }; 
(* Formulas: 
dplcl //InputForm 

(1 + vl4 + v23*vl4 - 2*vl4"2 - v23*vl4"2 - 2*v34 + 2*v23*v34 - 4*vl4*v34 + 
4*v23*vl4*v34 - 2*vl4"2*v34 + 2 *v2 3 *vl 4 " 2 *v34 + v34"2 + 2*v23*v34~2 - 
vl4*v34"2 - v23*vl4*v34"2 - v2 3* vl 4 " 2 * v3 4 " 2 + 
(1 + vl4) *Sqrt [vl4"2* (v23* (-1 + v34)~2 + 2*(1 + v34))'^2 + 

2*vl4*(2*(-l + v34)"2*(l + v34) + 2*v23~2*(-l + v34 ) " 2 * v3 4 * ( 1 + v34) + 
v23*(l - 18*v34"2 + v34"4)) + (1 + v34*(-2 + v34 + 2*v23*(l + v34))) 
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2])/{2*{l + v34)*{-l - (2 + v23)*vl4 + v34 + v23* (2 + vl4)*v34)) 
dplc2//InputForm 

(-1 + v23 - 2*vl4 + v23*vl4 + v23"2*vl4 + 2*v34 + 4*v23*v34 + 2*v23"2*v34 - 
2*vl4*v34 - 4*v23*vl4*v34 - 2 * v2 3 ~ 2 * vl 4 * v3 4 - v34"2 - v23*v34"2 + 
2*v23"2*v34~2 - v2 3* vl 4 * v3 4 ~ 2 + v2 3 " 2 *vl 4 *v34 " 2 + 
(1 + v23) *Sqrt [vl4"2* {v23* (-1 + v34)~2 + 2*{1 + v34))~2 + 

2*vl4*{2*{-l + v34)"2*{l + v34) + 2*v23"2*{-l + v34 ) " 2 * v3 4 * { 1 + v34) + 
v23*{l - 18*v34"2 + v34"4)) + (1 + v34*{-2 + v34 + 2*v23*{l + v34))) 
2])/{2*{l + v34)*{-l - (2 + v23)*vl4 + v34 + v23* (2 + vl4)*v34)) 
dplc3/ /InputForm 

{{-1 + v34)*(l + (2 + v23)*vl4 + v34 + v23*{2 + vl4)*v34) - 
Sqrt [vl4"2* {v23* (-1 + v34)"2 + 2*{1 + v34))"2 + 

2*vl4*{2*{-l + v34)"2*{l + v34) + 2*v23"2*{-l + v3 4 ) " 2 *v34 * { 1 + v34) + 
v23*{l - 18*v34"2 + v34"4)) + (1 + v34* (-2 + v34 + 2*v23*{l + v34)))'^ 
2])/{2*{-l - (2 + v23)*vl4 + v34 + v23* (2 + vl4)*v34)) 
dplc4// InputForm 
-1 

dp2cl// InputForm 

(1 + vl4 + v23*vl4 - 2*vl4"2 - v23*vl4"2 - 2*v34 + 2*v23*v34 - 4*vl4*v34 + 
4*v23*vl4*v34 - 2*vl4"2*v34 + 2 *v2 3 *vl 4 " 2 *v34 + v34"2 + 2*v23*v34~2 - 
vl4*v34"2 - v23*vl4*v34"2 - v2 3* vl 4 " 2 * v3 4 " 2 - 
(1 + vl4) *Sqrt [vl4"2* (v23* (-1 + v34)~2 + 2*{1 + v34))~2 + 

2*vl4*{2*(-l + v34)"2*(l + v34) + 2*v23"2*{-l + v34 ) " 2 * v3 4 * { 1 + v34) + 
v23*{l - 18*v34~2 + v34"4)) + (1 + v34*{-2 + v34 + 2*v23*{l + v34))) 
2])/(2*(l + v34)*(-l - (2 + v23)*vl4 + v34 + v23*{2 + vl4)*v34)) 
dp2c2// InputForm 

(-1 + v23 - 2*vl4 + v23*vl4 + v23"2*vl4 + 2*v34 + 4*v23*v34 + 2*v23"2*v34 - 
2*vl4*v34 - 4*v23*vl4*v34 - 2 * v2 3 ~ 2 * vl 4 * v3 4 - v34~2 - v23*v34"2 + 
2*v23"2*v34"2 - v2 3* vl 4 * v3 4 ~ 2 + v2 3 " 2 *vl 4 *v34 ~ 2 - 
(1 + v23) *Sqrt [vl4"2* {v23* (-1 + v34)"2 + 2*{1 + v34))"2 + 

2*vl4*{2*(-l + v34)"2*{l + v34) + 2*v23"2*{-l + v34 ) " 2 * v3 4 * { 1 + v34) + 
v23*{l - 18*v34~2 + v34"4)) + (1 + v34*{-2 + v34 + 2*v23*{l + v34))) 
2])/{2*{l + v34)*{-l - (2 + v23)*vl4 + v34 + v23* (2 + vl4)*v34)) 
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dp2c3/ /Input Form 

{{-1 + v34)*{l + (2 + v23)*vl4 + v34 + v23*{2 + vl4)*v34) + 
Sqrt [vl4"2* {v23* (-1 + v34)"2 + 2*{1 + v34))"2 + 

2*vl4*{2*{-l + v34)"2*{l + v34) + 2*v23"2*{-l + v3 4 ) " 2 *v34 * { 1 + v34) + 
v23*{l - 18*v34"2 + v34"4)) + (1 + v34* (-2 + v34 + 2*v23*{l + v34))) 
2])/{2*{-l - (2 + v23)*vl4 + v34 + v23* (2 + vl4)*v34)) 
dp2c4// Input Form 
-1 

* And 11*12 == 
lltl2//InputForm 

-{{-1 + ph)*(l + v34)"2*(-l - (2 + v23)*vl4 + v34 + v23*{2 + vl4)*v34))/ 
{4*{-l + v23*vl4) "2* (-1 + v34)*v34) 



